Optimal Angular Resolution for Face-Symmetric Drawings 



David Eppstein and Kevin A. Wortman 
Department of Computer Science 
Universitiy of California, Irvine 

{eppstein, kwortman}@ics. uci.edu 

July 31, 2009 

O 
O 

<N 

Abstract 

l—j Let G be a graph that may be drawn in the plane in such a way that all internal faces are 

I centrally symmetric convex polygons. We show how to find a drawing of this type that max- 

imizes the angular resolution of the drawing, the minimum angle between any two incident 
edges, in polynomial time, by reducing the problem to one of finding parametric shortest paths 
in an auxiliary graph. The running time is at most 0(t 3 ), where t is a parameter of the input 
graph that is at most 0(n) but is more typically proportional to y/n. 

,°, 1 Introduction 

Angular resolution, the minimum angle between any two edges at the same vertex, has been rec- 
^j- ognized as an important aesthetic criterion in graph drawing since its introduction by Malitz and 

Papakostas in 1992 [10]. Much past work on angular resolution in graph drawing has focused on 
bounding the resolution by some function of the vertex degree or other related quantities, rather 
than on exact optimization; however, recently, Eppstein and Carlson [3] showed that, when draw- 
ing trees in such a way that all faces form (infinite) convex polygons, the optimal angular resolu- 
Q\ tion may be found by a simple linear time algorithm. The resulting drawings have the convenient 

property that the lengths of the tree edges may be chosen arbitrarily, keeping fixed the angles se- 
lected by the optimization algorithm, and no crossing can occur; therefore, one may choose edge 
lengths either to achieve other aesthetic goals such as good vertex spacing or to convey additional 
H information about the tree. 

In this paper, we consider similar problems of calculating edge angles with optimal angular 
resolution for a more complicated class of graph drawings: planar straight line drawings in which 
each internal face of the drawing is a centrally-symmetric convex polygon. In previous work [6], 
we investigated drawings of this type, which we called face-symmetric drawings. We characterized 
them as the duals of weak pseudoline arrangements in the plane, and described an algorithm 
that finds a face-symmetric drawing (if one exists) in linear time, based on an SPQR-tree decom- 
position of the graph. Graphs with face-symmetric drawings are automatically partial cubes (or, 
equivalently, media [7]), graphs in which the vertices may be labeled by bitvectors in such a way 
that graph distance equals Hamming distance; see [7] for the many applications of graphs of this 
type. 
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Figure 1: A hyperbolic line arrangement (left), and its dual squaregraph (right), drawn by our 
previous un-optimized algorithm with all faces rhombi. 

Figure[T]depicts an example, a hyperbolic line arrangement with no three mutually intersecting 
lines, the intersection graph of which requires five colors, as constructed by Ageev [1], and the 
planar dual graph of the arrangement. The dual of a hyperbolic line arrangement with no three 
mutually intersecting lines, such as this one, is a squaregraph, a type of planar median graph in 
which each internal face is a quadrilateral and each internal vertex is surrounded by four or more 
faces [2,5]. Any squaregraph may be drawn in such a way that its faces are all rhombi, by our 
previous algorithm [6], and the drawing produced in this way is shown in the figure. However, 
note that some rhombi have such sharp angles that they appear only as line segments in the figure, 
making them difficult to view. Other rhombi, even those near the edge of the figure, are drawn 
with overly wide angles, making them very legible but detracting from the legibility of other 
nearby rhombi. Thus, we are led to the problem of spreading out the angles more uniformly across 
the drawing, in such a way as to optimize its angular resolution, while preserving the property 
that all faces are rhombi. 

As we show, this problem of optimizing the angular resolution of a face-symmetric drawing 
may be solved in polynomial time, by translating it into a problem of finding parametric shortest 
paths in an auxiliary network. In the parametric shortest path problem, each edge of a network 
is given a length that is a linear function of a parameter A; substituting different values of A into 
these functions gives different real weights on the edges and therefore different shortest path 
problems [8]. The number of different shortest paths formed in this way, as A varies, may be 
superpolynomial [4], but fortunately for our problem we do not need to find all shortest paths 
for all values of A. Rather, the optimal angular resolution we seek may be determined as the 
maximum value of A for which the associated network has no negative cycles, and the drawing 
itself may be constructed using distances from the source in this network at the critical value of 
A. The linear functions of our auxiliary network have a special structure — each is either constant 
or a constant minus A — that allows us to apply an algorithm of Karp and Orlin [9] for solving this 
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Figure 2: Converting a face-symmetric drawing into a weak pseudoline arrangement, by drawing 
line segments connecting opposite pairs of edge midpoints within each face. From [6]. 

variant of the parametric shortest path problem, and thereby to optimize in polynomial time the 
angular resolution of our drawings. 

2 Drawings with symmetric faces 

The algorithm of [6] can be used to find a non-optimal face-symmetric drawing, when one exists, 
in linear time; therefore, in the algorithms here we will assume that such a drawing has already 
been given. We summarize here the relevant properties of the drawings produced in this way. 

In a face-symmetric drawing, any two opposite edges of any face must be parallel and of equal 
length. The transitive closure of this relation of being opposite on the same face partitions the 
edges of the drawing into equivalence classes, which we call zones; any zone consists of a set of 
edges that all have equal angles and equal lengths. (Note, however, that edges in different zones 
may also have equal angles and equal lengths.) The collection of line segments connecting oppo- 
site pairs of edge midpoints within each face forms a collection of curves which can be extended 
to infinity to form a weak pseudoline arrangement; each zone is formed by the drawing edges that 
cross one of the curves of this arrangement. This construction is depicted in Figure [2] 

Thus, the positions of all the vertices of the drawing are determined, up to translation of the 
whole drawing, by a choice of a vector for each zone that specifies the direction and length of the 
zone's edges. For, if an arbitrarily chosen base vertex has its placement fixed at the origin, the po- 
sition of any other vertex v must be the sum of the vectors corresponding to the weak pseudolines 
that separate v from the base vertex. All vertex positions may be computed by performing a depth 
first search of the graph, setting the position of each newly visited vertex v to be the position of v's 
parent in the depth first search tree plus the vector for the zone containing the edge connecting v 
to its parent. 

The algorithm from [6] determines the vectors for each zone as follows. Associate a unit vector 
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with each end of each of the pseudolines dual to the drawing, in such a way that these unit vectors 
are equally spaced around the unit circle in the same cyclic order as the order in which the pseu- 
doline ends extend to infinity. The vector associated with each zone is then simply the difference 
of the two unit vectors associated with the corresponding pseudoline's ends, normalized so that 
it is itself a unit vector. 

As shown in [6], this choice of a vector for each zone leads to a face-symmetric drawing with- 
out crossings. Additionally, it has two more properties, both of which are important and will 
be preserved by our optimization algorithm. First, if the planar embedding chosen for the given 
graph has any symmetries, they will be reflected in the dual pseudoline arrangement, in the choice 
of zone vectors, and therefore in the resulting drawing. 

Second, although it may not be possible to draw the given graph in such a way that the outer 
face is convex, its concavities are all mild. This can be measured by defining the winding number of 
any point p on the boundary of the drawing, with respect to any other point q also on the boundary, 
as the sum of the turning angles between consecutive edges along a path counterclockwise around 
the boundary from q to p. For any simple polygonal boundary, the winding number from p to q 
is 2-7T minus the turning angle from q to p. In a convex polygon, all winding numbers would lie in 
the range [0, 2n]; in the drawings produced by this method, they instead lie in the range [— n, 3tt]. 
That is, intuitively, the sides of any concavity may be parallel but may not turn back towards each 
other. It follows from this property that the vectors of each zone may be scaled independently of 
each other, preserving only their relative angles, and the resulting drawing will remain planar [6]. 
In particular, the step in the algorithm of [6] in which the zone vectors are normalized to unit 
length leads to a planar drawing. A similar constraint on winding numbers was used in the 
algorithm of [3] for finding the optimal angular resolution of a tree drawing, and again led to the 
ability to adjust edge lengths arbitrarily while preserving the planarity of the drawing. Indeed, 
tree drawings may be seen as a very special case of face-symmetric drawings in which there are 
no internal faces to be symmetric. 

Thus, we may formalize the problem to be solved, as follows: given a face-symmetric drawing 
of a graph G, described as a partition of the edges of G into zones and a zone vector for each zone, 
we wish to find a new set of zone vectors to use for a new drawing of G, preserving the relative 
orientation of any two edges that meet at a vertex of G, maintaining the constraint that all winding 
numbers lie in the range [— ir, 3ir], and maximizing the angular resolution of the resulting drawing. 

3 Algorithmic results 

In this section we describe a polynomial time algorithm for the problem at hand. As described 
above, the core of our algorithm is a routine that maps a planar face-symmetric drawing G to an 
auxiliary graph A, such that the auxiliary graph may be used as input to a parametric shortest 
paths algorithm of Karp and Orlin [9], and the output of that algorithm describes a drawing G' 
isomorphic to G and with maximal angular resolution. The algorithm of [9] runs in polynomial 
time, so what remains to be seen is the correctness and running time of our mapping routine. 

Table [l] summarizes the intuitive relationships between concepts in the drawings G and G' and 
their representation in the auxiliary graph A. As shown in the table, edge zones correspond to 
auxiliary graph vertices; angular resolution corresponds to the parametric variable A; the change 
in angle of zones between G and G' corresponds to lengths of shortest paths in A; and constraints 
on the angles are represented by edges in A. Figure [3] shows a small example input graph G with 
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Concept in drawing 


Representation in auxiliary graph 


Zone Zi 


Vertex Vi 


Angular resolution 


Parametric variable A 


a is feasible angular resolution 


No negative cycles when A < a 


Difference between angles of vectors for zone 
Zi in unoptimized and optimized drawings 


Shortest path distance d(v.i) from s to v- L 


Angle between z\ and Zj is > A 


Edge with weight Og(z{) — 9c{zj) — A 


Interior faces are convex 


Edge with weight it + 9c(zi) — Oc( z j) 


Exterior boundary is mildly convex 


Opposing edges with weight w(vi,Vj) = 37r + 
QG(zi) ~ Sc(zj) and w{vj,Vi) = tt + 9 G {zj) - 
G {zi) 



Table 1: Relationships between drawing concepts and auxiliary graph components. 




Figure 3: Example input graph G, using degrees rather than radians. 

five zones, and Figure |4] shows the corresponding auxiliary graph A. 

We now define our notation formally and show that our reduction is correct. 

Definition 3.1. Let 

• Gbe the planar face-symmetric graph and its drawing given as input; 

• Z = {zq, z\, . . .} be the zones ofG; 

• t = \Z\ be the number of zones in G; 

• Abe a weighted, directed, auxiliary graph, such that every zone zi G Z corresponds to a vertex Vi in 
A, and every edge of a has weight w(vi, Vj) = b V)W — m V)W ■ X, where 5 G 1 and m v>w £ {0, 1} are 
per-edge constants, and A is a graph-wide variable; 

• sbea special start vertex in A; 

• A have edges (vi, Vj) and (vj,Vi)for every pair of vertices Vi,Vj; 
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Figure 4: Auxiliary graph A for the input shown in Figure|3j again using degrees. Only the lowest- 
weight edge between any pair of vertices is shown. 

• A* be the largest value of X such that A contains no negative cycles; 

• d{vi) be the weight of the shortest path from s to v% in A when A = A*; 

• 9 G (zi) be the angle assigned to edges of zone z% in G, in radians counterclockwise from the x-axis; 

• Z G (i, j) = Qc(zj) — Oc(zi) be the counterclockwise turning angle from 9a{zi) to 9 G {zj) in G; 

• for every pair of distinct zones Zi, zj, such that 9 G {zi) < 9c(zj) and there exists a face in G with 
some edge from Zi and some edge from Zj, A contains the following edges: 

- an edge (vj,Vi) with weight w(vj,Vi) = Og(zj) — 9a{zi) — K 

- if a corresponding angle in G is interior, an edge (v{,Vj) with weights w(vi, Vj) = tt + Oc(zi) — 
8a(zj), and 

- if a corresponding angle in G is exterior, opposing edges with weights w(vi,Vj) = 3tt+9g{zi) — 
9c{zj) andw(vj,Vi) = it + 9 G {zf) - 9 G {zi); 

• G' be the output drawing of G; 

• 9qi{z{) = d(vi) + 9c(zi) be the angle assigned to edges in zone Z{ in the output drawing G' , in 
radians counterclockwise from the x-axis; 

• and Z G i(i,j) = 9 G i{zf) — 9 G >{zi) be theangle between Z{ and Zj in G', analogous to Z G (i,j) in G. 
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Lemma 3.2. If A contains any edge from v,- t to Vj with weight w(vi, vf) = W, then d(vj) < d(vi) + W. 

Proof. The shortest path to V{ in A, followed by the edge (vi, vf), is a path to Vj with total weight 
d(vi) + W. Since d is defined for X = X*, A has no negative cycle; so the shortest path to Vj has 
weight d(vj) < d(vi) + W. □ 

Lemma 3.3. G' has angular resolution A*, and every interior face of G' is non-concave. 

Proof. Let e% and ej be any pair of edges in G that form an angle on some interior or exterior face. 
If e« has the lesser absolute angle measure, then by construction there exists an edge (vj, vi) with 
weight w(vj,Vi) = 9 G (zj) — 9 G (zi) — A*, where V{ and vj are the zones for a and ej, respectively. 



Thus by Lemma 3.2 





d{vi) 


< 


d(vj) + w(vj,Vi) 




d(vi 


)-d(Vj) 


< 


(Oaizj) - e G { Zi ) 


- A*) 


-d(vi 


) + d{Vj) 


> 


-o G {zj) + e G ( Zi ) 


+ A* 


zj) - d(vi) 


- G {Zi) 


> 


X* 








> 


X* 








> 


A*; 





so every corner angle in G' has measure no less than A* . 

Now suppose Zeiej is interior. Then by construction there also exists an opposing edge (v{, Vj) 



with weight w(vi, vj) = ir + G (zi) — 9 G (zj). So again by Lemma 3.2 



d(vj) < d( Vi ) + (tt + 9 G ( Zi ) - e G ( Zj )) 
d(vj) + 9 G (zj) - d{vi) -9 G {zi) < tt 
(6 GI (z :i ))-(9 G/ (z l )) < TT 

which implies that G' contains no concave face. □ 
Lemma 3.4. The winding number of any p and q on the boundary of G' is in the range [—tt, 3tt]. 



Proof. As in Lemma |33} let and ej be edges in G such that 9 G (i) < 9 G (j). If ej and ej form an 
angle on the boundary of G, then the corresponding vertices Vi and Vj in A have opposing edges 
with weights w(vi,Vj) = 37r + 9 G (zi) — 9 G (zj) and w(vj, vi) = tt + 9 G (zj) — 6 G (zi). Then by algebraic 



manipulations symmetric to those in Lemma 3.3 



-TT < Z.G'(i,j) < 37T. 

□ 

Theorem 3.5. The output graph G' is a planar face symmetric drawing isomorphic to the input drawing 
G, such that all boundary winding numbers lie in the range [—tt, 3tt], and the angular resolution of G' is 
maximal for any such drawing. Given G, G' may be generated in 0(t 3 ) time. 
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Proof. Identifying the zones of G and constructing A may be done naively in 0(t 2 ) time. The 
value A* and corresponding distances d(vi) for every vertex v% in A may be reported by invoking 
the algorithm of Karp and Orlin. This algorithm runs in 0(n 3 ) time, where n is the number of 
vertices in the input; our graph A has t vertices, so this step takes 0(t 3 ) time. The d(vi) distances 
define the zone angles 9q', allowing the drawing G' to be output in linear time as described in 
Section |2j 

Our method of choosing new angles for the zone vectors implies that G' is isomorphic to G 
and every interior face of G' is centrally symmetric. By Lemmas |3.3| and 3.4 every interior face of 
G' is convex, and the exterior boundary of G' satisfies the convexity constraint. 

Finally we must show that G' has maximal angular resolution. Let H be any correct drawing 
of the graph G, and let H have angular resolution A#. Then A has no negative cycles for A = Xh- 
for, suppose we replace every weight w(v{, vj) with w'(vi, vj) = w(vi,Vj) + 9h{vi) — @H(vj), where 
9ii{vk) is the angle assigned to zone zu in H. Any cycle contains an equal number of +9n(vk) and 
—&H(vi) terms for each Vi in the cycle, so our transformation does not change total cycle weights, 
and hence preserves negative cycles. Each edge has nonnegative weight, so no negative cycles 
exist in A. We use the Karp-Orlin algorithm to find the largest value A* for which A contains 
no negative cycle, and by Lemma |3.3| the output graph G' has angular resolution A* . Thus the 
angular resolution of G' is greater than or equal to that of any correct drawing H. □ 



4 Experimental results 

We implemented a simplified version of the algorithm described in Section [3] in the Python lan- 
guage. Our simpler algorithm performs a numerical binary search for A* rather than an imple- 
mentation of the Karp-Orlin parametric search algorithm. We make binary search decisions by 
generating the auxiliary graph as described in Section |3j substituting in the value of A to obtain 
real-valued edge weights, then checking for negative cycles with a conventional Bellman-Ford 
shortest paths computation. This algorithm runs in 0(t s log W) time, where t is the number of 
pseudolines (zones) in the drawing and W is the number of bits of numerical precision used. Wall 
clock time ranged from a matter of seconds for small t to roughly two minutes for t = 220 and 
W = 64. 

Figure [5] shows a sample of input and output for our implementation. The input is a correct, 
though suboptimal, planar face-symmetric drawing of a graph with 15 pseudolines generated by 
the algorithm presented in [6]. The angular resolution of the output is visibly improved. Figure 
[6] shows the output of the algorithm when the 220-pseudoline graph shown in Figure [T] is used 
as input. That unoptimized drawing has angular resolution 2n/U0 ~ .0571 radians. The left 
output drawing was produced by our optimization algorithm using all correctness constraints, 
and has angular resolution 2it/75 as .0838 radians. The right drawing is the result of removing the 
constraint that concavities have opening angles of at least tt radians, which may in general result 
in a nonplanar drawing, but in this case yields a planar drawing with an even greater angular 
resolution of 27r/30 rj .209 radians. 
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Figure 5: Example input (left) and output (right) with 15 pseudolines. 



5 Conclusion 

We have described two algorithms for generating face-symmetric drawings with optimal angular 
resolution. The first algorithm runs in strictly cubic time and uses a subsidiary parametric short- 
est paths algorithm as a black-box subroutine. The second algorithm runs in pseudopolynomial 
time and relies on the simpler Bellman-Ford shortest paths algorithm. We implemented the latter 
algorithm and found that it generates output that is both numerically and visibly improved over 
unoptimized drawings. 

Finally we offer the following possible directions for future research: 

• We choose the edge length for each zone arbitrarily; can choosing them more carefully lead 
to improved legibility? 

• Can this approach be applied to related types of drawings, such as the projections of high 
dimensional grid embeddings also studied in [6]? 

• Our algorithm respects a fixed embedding. What about allowing for different embeddings, 
e.g. flipping parts of the graph at articulation vertices? Is it still possible to optimize angular 
resolution efficiently in this more general problem? 

• What about other angle optimization criteria, such as those defined by all angles in the draw- 
ing, rather than merely the sharpest angle? 

References 

[1] A. A. Ageev. A triangle-free circle graph with chromatic number 5. Discrete Mathematics, 
152:295-298,1996. 

[2] H.-J. Bandelt, V. Chepoi, and D. Eppstein. Combinatorics and geometry of finite and infinite 
squaregraphs. In preparation, 2009. 



9 



Figure 6: Output of our implementation for the graph shown in Figure [TJ safely optimized (left) 
and unsafely optimized (right). 

[3] J. Carlson and D. Eppstein. Trees with convex faces and optimal angles. In M. Kaufmann 
and D. Wagner, editors, Proc. 14th Int. Symp. Graph Drawing, volume 4372 of Lecture Notes in 
Computer Science, pages 77-88. Springer- Verlag, 2006. 

[4] P. D. Carstensen. Parametric cost shortest path problems. Unpublished Bellcore memo, 1984. 

[5] V. Chepoi, F. Dragan, and Y. Vaxes. Center and diameter problem in planar quadrangulations 
and triangulations. In Proc. 13th ACM-SIAM Symp. Discrete Algorithms (SODA), pages 346- 
355,2002. 

[6] D. Eppstein. Algorithms for drawing media. In J. Pach, editor, Proc. 12th Int. Symp. Graph 
Drawing (GD 2004), volume 3383 of Lecture Notes in Computer Science, pages 173-183. Springer- 
Verlag, 2004. 

[7] D. Eppstein, J.-C. Falmagne, and S. Ovchinnikov. Media Theory. Springer- Verlag, 2007. 

[8] D. Gusfield. Parametric combinatorial computing and a problem of program module distri- 
bution. /. ACM, 30(3)551-563, 1983. 

[9] R. M. Karp and J. B. Orlin. Parametric shortest path algorithms with an application to cyclic 
staffing. Discrete Appl. Math., 3(l):37-45, 1981. 

[10] S. Malitz and A. Papakostas. On the angular resolution of planar graphs. In Proc. 24th ACM 
Symp. Theory of Computing, pages 527-538, 1992. 



10 



